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Abstract 

A recently developed quark model approach to pseudoscalar meson photo- 
production is extended to electroproduction process for the r] meson in the 
kinematics of momentum transfer < 4 (GeV/c)^ and total center of mass 
energy W < l.Q GeV. Existing data are well reproduced and the roles of the 
5ii(1535) and Di3(1520) resonances are closely investigated. In the study of 
the longitudinal excitation of the 5ii(1535) resonance, a reliable constraint 
on the 5*11(1535) properties is obtained by cleanly removing the electromag- 
netic transition from the 7(t,)P S'ii(1535) rjp amplitude. Thus, the fitted 
quantities can be determined with an uncertainty of about 15%. This could 
be the first direct constraint on the 5ii(1535) properties in theory. 
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I. INTRODUCTION 



Our knowledge about the internal structures of baryon resonances, which generally belong 
to non-perturbative QCD phenomena, are still far from complete. Nowadays, with the 
availabilities of high intensity photon and electron beams at JLab, ESRF, MAMI, ELSA, 
and SPring-8, the baryon resonances thus can be systematically investigated via meson 
photo- and electroproduction. This initiates various theoretical efforts, through which one 
expects that our knowledge about these non-perturbative phenomena can be established on 
a more solid and fundamental ground. 

In the late 1990 's, experimental data with unprecedent accuracies were available for the rj 
meson photoproduction which allowed a close study of resonance excitations in 7(i,)p rjp. 
For example, the photoproduction data for the cross sections and single polarization ob- 
servables obtained at MAMI |], ELSA § and ESRF (GRAAL) gi, had allowed to es- 
tablish that the rj photoproduction reaction mechanism was dominated by the S'ii(1535) 
excitation, while the i5i3(1520) played a small but non-negligible role ||5|-[rT|. The measure- 
ments for polarization observables also made it possible to single out the contribution of 
the Fi5(1680) resonance However, detailed properties concerning the rjNN* couplings 
and even the rjNN coupling still could not be well-constrained at that time, including the 
total decay width of the dominant S'ii(1535), and its branching ratio into riN. Such a 
situation somehow stalked theoretical progresses due to lack of constrained inputs for the 
modellings, although various QCD inspired phenomenologies, such as the quark model ap- 
proaches PJT^-[l4|,^, and formalisms based on the meson- nucleon degrees of freedom using 
isobaric descriptions ||7|j9|JTl|l , or multipole analyses P,|lO|, were proposed to investigate the 
roles played by the intermediate resonances in the rj photoproduction. 

The situation was changed significantly with the electroproduction data available from 
JLab [|l^, which brought rich complementary information about the 5*11(1535) properties 
and the D-wave influences, due to the presence of the longitudinal photon excitations. Such 
an observable is not only valuable for model-selection, but also useful for disentangling 
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model-dependent or model-independent aspects within a phenomenology. 

In this paper, we shall study the rj meson electroproduction in a constituent quark model 
with a chiral effective Lagrangian for the quark-meson [r]) couplings, and focus our attentions 
on the dominant 5*11(1535) and Di^{1520) resonance. The formalism is discussed in Sec. II. 
In Sec. Ill, we present detailed results produced by our model fits. This will allow us to 
proceed to a successful constraint on the 5*11(1535) properties. The conclusions are given in 
Sec. IV. 



II. THEORETICAL FRAME 

In this work, we follow the scheme of a recently developed quark model approach to the 
pseudoscalar meson photoproduction [|T2|, and extend it to the rj meson electroproduction 
in the region of the cm. energy W = 1.535 GeV. Compared to an isobaric approach, in this 
model the intermediate baryon resonances can be systematically included in the formalism. 
Starting from the NRCQM [|16|, only a limited number of adjustable parameters will appear 



in the model. 



A. Formalism 



For the meson interaction vertex, a QCD inspired effective Lagrangian [|T7l is introduced 
to account for the quark-meson interaction: 

where ipj (ipj) is the jth quark (anti-quark) field in the nucleon, and (prj represents the rj meson 
field. At quark level, the transition matrix element Ai ji can be expressed as the sum over the 
t-, s- and w-channel transitions, i.e. Aiji = M.^fi + A^jj + A^jj. Since M.^fi is proportional to 
the charge of the outgoing meson, it vanishes in the neutral meson production. In the s- and 
M-channel, contributions from a complete set of intermediate baryon resonances are included 
in the quark model SU(6)C>?'0(3) symmetry limit. It can be seen easily that due to the 



isospin conservation, isospin 3/2 states cannot contribute in this channel. Nevertheless, the 
Moorhouse selection rule eliminates the states of representation [70, '^8] from contributing 



to the proton target reaction [|T^. Thus, eight isospin-1/2 resonances corresponding to the 
harmonic oscillator quantum numbers n=l and 2, and the nucleon pole terms, are to be 
included explicitly. States with n > 2 are treated as degenerate in the quantum number n. 

Adopting Lorentz gauge = 0, the longitudinal virtual photon polarization vector 

is defined as 

= -^(|k|, 0, 0, u,) , (2) 

where |k| and uj^ are the momentum and energy of the virtual photon in the meson-nucleon 
cm. frame, respectively. Gauge invariance requires that the longitudinal transition operator 
is proportional to v^Q^/|k|. Therefore, the longitudinal transition vanishes in the real photon 
limit. The longitudinal electromagnetic interaction is defined as 

Hem = ^0^0 - £^3-^3 , (3) 

where the electromagnetic current components Jq and J3 for a three-quark system can be 
expressed as 

Jo = E e,e^-^ , (4) 

and 

Ja = E e,a,e^-^ , (5) 

with Cj denotes the charge operator of the jth quark, and aj, the Dirac matrix. Given the 
Hamiltonian of the three-quark system described by only two-body interactions |T^, the 
current conservation gives 

{Nfij^m = ^^^{NfiMm , (6) 

where ks = |k| and the energy conservation gives Ej — Ei = uo^. The longitudinal interaction 
can then be written as 



H 



(7) 



Substituting the into the above equation, one obtains the gauge-invariant relation: 



j=l I I 



(8) 



We adopt the following nonrelativistic expansion of Jq for the longitudinal transition 

1 



J, 



NR 



-E 



i a J 



4Mt m,- 



mi ■ 



ejk X piC*'"''"^ - e^k x Pje"" ""' } 

j<i ^^"J ""3 

where the first term describes the cm. motion of the three-quark system and is commonly 
used in the longitudinal transitions. The second and third terms are spin-orbit and non- 
additive terms, which are believed to be the leading order relativistic corrections to the 
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transition operators We shall neglect the spin-orbit and non-additive terms in this 

work. It can be seen later that this approximation will result in an over-estimation of 
the 5*11 (1535) longitudinal excitation cross section. Empirically, a coefficient for the Sn 
longitudinal transition amplitude will be introduced. 

The transverse transition amplitude for the photoproduction have been derived in 



Ref. [12|, in which the transition amplitudes for the S'ii(1535) and i5i3(1520) are 
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where /ig is the reduced mass of two quarks and equals mg/2 in the t] production for m^ = 
m^ = mq. The overall coupling is related to the rj meson radiative decay constant by 



assuming the validity of the Goldberger-Treiman relation p2[ for the rjNN couplings: 
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with = g'^j^j^/iTT, where = 1 is given by the NRCQM in the SU(6)®0(3) symmetry 
hmit. 

Similarly, the amplitudes for other resonances in the quark model can be derived. In 
this paper, we concentrate on the kinematics from threshold to W ~ 1.54 GeV, where 
the available experimental data allow a close study of the 511(1535) and Di3(1520). One 
of the advantages is that one can neglect the SU(6)®0(3) symmetry violations for other 
excited states. Therefore, we could only introduce an SU(6)(8>0(3) beaking coefficient for 
the 5*11(1535) and Di3(1520). In another word, except for the 5ii(1535) and Di3(1520), the 
relative strengths among other resonances will be constrained by the quark model. In this 
way, we avoid to introduce too many parameters in the model. This treatment is different 



from that in Refs. P and Another feature in this calculation is that the quark model 
parameters are also taken into account. The harmonic oscillator strength a is treated as 
a free parameter, of which a value within the range of quark model validity could be an 
essential test of the self-consistence of this model. Moreover, the 5ii(1535) total width F^^^ 
and its partial decay branching ratio 6,, into riN will be treated as free parameters as well. 

B. Kinematics and observables 

The kinematics of the meson exclusive electroproduction has been discussed in the lit- 
erature [24,^1 . Here, we directly relate the kinematics to our convention of the transition 



amplitude. Generally, the cross section of the meson electroproduction can be written as 
follows: 

. (13) 



dE'jn2dn* dn* 

where, da/dQ* represents the cross section for the virtual-photon- nucleon scattering in the 
meson-nucleon cm. frame, while the contribution from the lepton current can be factorized 
into a factor r^, which is known as the virtual photon flux 
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where is the electromagnetic fine structure constant, Eg and E'^ are the energies of the 
initial and scattered electrons, respectively, in the lab system, and Ke = (s — M^)/2M is 
the equivalent energy of the virtual photon as a real photon in the lab system. 
The virtual photon polarization parameter, e, is defined as 

£ = (l + 2tan2(|)|ko|VQ')-'- (15) 

where ko is the virtual photon momentum in the lab system. 

The cross section for the meson production by the virtual photon in the meson-nucleon 
cm. frame can be expressed as 

dn^^dn^ + 'dn^-'ln^''' ^"^^ ' + '^In^''' ' ^^^^ 

where 0* is the azimuthal angle between the r]N scattering plane and the (e, e') scattering 
plane, and 
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In the above equations, ^ represents the phase space factor in the virtual-photon-nucleon 
excitations, and has the following expression: 

ae Mn |q| 
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where |q| = (s — Mj^ — m^Y — AM^m^ ^ /(21^) is the meson momentum in the meson- 
nucleon cm. frame. The amplitude, H'^^\ is defined as 

n''' = E KxAxJ > (22) 
A1A2 



where ^a^a2 ^'^^ transition amplitudes for the meson production via the virtual photon, 
A = 0, ±1 denotes the spin polarizations of the incoming virtual photon and, Ai = ±| 
and A2 = ±^ for the initial and final state nucleons, respectively. In the pseudoscalar 
meson electroproduction, there are two independent amplitudes in the longitudinal photon 
transitions, and four independent amplitudes in the transverse photon transitions. The latter 
ones can be related to the traditional CGLN amplitudes or the helicity amplitudes ^6 . 



When the calculations are extended to the electroproduction, it becomes highly relativis- 
tic, and a Lorentz boost factor must be introduced into the spatial integrals. We follow the 



prescription of Foster and Hughes p3 boost momentum for the spatial integrals in an 
equal velocity frame (EVF), which is very close to the Breit frame. In the EVP, the Lorentz 
boost factor for the virtual photon interaction is defined as 

where the relation of the momentum k with the in the EVF is 



,^(WF)=(^^l^^+«!(15i±^. (24) 



With the boost factor, the spatial integral concerning the photon excitation is boosted as 

R{k) \r{—). (25) 

Ik Ik 

The same prescription can be adopted for the resonances decaying into the final state r\ 
meson and proton: 

7, = fl + TTiT^Z^V ' (26) 



where 



Therefore, the spatial integrals in the transition matrix elements from the initial states 
to the final states can be expressed as 
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Compared to the previous study of the t] meson photoproduction , in which the Lorentz 
boost factors, 7^ = Ei/M^- and 7^ = Ef/Mf, in the cm. frame of the meson-nucleon system 
were adopted, the Lorentz boost factors in the EVF turn out to be more successful in 
reproducing the data. The failure of the Lorentz boost factor 7^ in the meson-nucleon cm. 
frame can be seen more clearly by rewriting it in the electroproduction: 



Since 4sM^/(s + M^)^ < 1, the rest part plays a role as a dipole in the nucleon's elec- 
tromagnetic form factor. Note that, the cm. energy W = y/s is fixed at M5^^, therefore, 
(s + M^) is a constant. In contrast to the Lorentz boost factor in the EVF, an ad hoc dipole, 
1/(1 + (5^/0.71)^, was widely used for the nucleon's electromagnetic form factor in the litera- 
ture ||2^. Such a Q^-dependence results in a fast drop-down behavior, which underestimates 
the data significantly at high regions. 

To show the difference between the two Lorentz boost schemes, we present the Q^- 
dependence of the I/7I factor calculated at W = 1.535 GeV by Eqs. (solid curve) and 
p9| (dashed curve) in Fig. |l|. It shows that in the EVF, electromagnetic transition is boosted 
rather moderately than in the meson-nucleon cm. frame. The difference between these two 
frames reveals non-covariance of this approach. As the boost goes large, a more realistic 
prescription is needed. 

For the meson interaction vertex, the boost factor at W = 1.535 GeV in both frames 
is a constant, and possesses very close values, i.e., 1/7^ = 0.99 in the EVF and 0.96 in the 
meson-nucleon cm. frame. Such a feature is due to the quite heavy mass of the 77 meson 
and its relatively small on-mass-shell momentum. 




(29) 



III. RESULTS AND DISCUSSIONS 



In this Section, we present numerical results for various photo- and electroproduction 
observables. Parameters are fitted by experimental data. Meanwhile, an analytical relation 
is deduced in the study of the 5*11(1535) excitations. 

A. Database and adjustable parameters 

The newly published data cover the kinematics 1.490 <W< 1.615 GeV at =2.4 



and 3.6 (GeV/c)^, have greatly improved the experimental status of the t] meson photo- 
and electroproduction. Along with three other datum sets: old electroproduction data at 



W= 1.535 GeV with =0.20, 0.28 and 0.40 (GeV/c)^ |9[, and recent photoproduction 
cross section data from Mainz |]l| and GRAAL |Q from rjN threshold to i?^"^ < 0.9 GeV 
(corresponding toW = 1.6 GeV), we obtain a large database covering the complete kinemat- 
ics of the 5*11(1535) and i5i3(1520) excitations. This opportunity will allow us to concentrate 
on these two resonances, whose couplings to the r]N channel have not been well known. Note 
that, the Mainz data cover the energy range from threshold to ^^'^''=0.789 GeV, while the 
GRAAL data go to higher energies: 0.714 < E;^"'' < 1.1 GeV. For GRAAL data we use a 
subset restricted to El^""^ < 0.9 GeV. 

The complete database contains 677 differential cross-sections and has been used to fix 
7 adjustable parameters in this approach. The results, obtained using the CERN MINUIT 
minimization code, are given in Table |. Below, we summarize the main features of each 
parameters based on the numerical investigations. 

The TjNN coupling constant (a^ = g'^^^/An) is not a well- determined quantity in the 
literatures. Here, is an overall factor and the numerical fits give = 0.10 ±0.01. Recent 
studies |]8| j30|j31[] suggest that g^j^j^/in is much smaller than unit. 

The commonly used value for the second free parameter, namely, the harmonic oscillator 
strength a varies between 330 and 430 MeV in the literatures. In the numerical fitting. 
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a rather moderate value, a = 384.5 MeV, is extracted. With the quark model spatial 
wavefunctions boosted in the EVF, Such a value avoids the introduction of an ad hoc dipole 
or monopole in the form factor 0. Then, the exponential factor from the spatial integrals, 
g-(k^+q^)/6a2^ wiU play a role like a form factor after being boosted. 

The 5*11 (1535) 's total width F^^^ and its t]N decay branching ratio 6^ have been very 
important in the study of its nature. However, both experimental values and theoretical pre- 
dictions for these two quantities have been very contradictory in the literatures ||T1, [T5| , |32| , |33 



The f] meson photo- and electroproduction at the Su mass region show great sensitivities to 
these two quantities. Much stronger constraints from the electroproduction data are found in 
this calculation. Treating the two quantities as free parameters, we extract r5^-^ = 143.3±0.2 
MeV and 6r;('S'ii) = 0.55 ± 0.01. These results are very close to the most recent extraction 



from the JLab data . For other resonances, their total widths and partial decay branch- 
ing ratios are found not to be sensitive in this kinematical region. Therefore, we use the 
total widths from the Particle Data Group (PDG) in the calculation. The partial decay 
channels for those resonances are restricted to the vrA^ and t]N, with branching ratios fixed 
at 6^=90% and 6^=10%. As a sensitivity test, we also use bT,=99% and 6^=1% and found 
that the numerical results do not show significant dependence on these branching ratios due 
to the dominance of the S'ii(1535). 

The fifth and sixth parameter, Cs^ and C^i^g, are strength coefficients introduced for the 
5*11(1535) and -Di3(1520) amplitudes to take care of the breaking of the SU(6)®0(3) symme- 
try. As discussed in Sec. II, eight resonances with harmonic oscillator shell n < 2 are to be 



^ In Ref. |27|, the authors adopted a larger value for the harmonic oscillator strength, a = 458.2 



MeV. Meanwhile, a monopole form factor, Fg{Q'^) = l/(l + (5^/0.8), had to be adopted to produce 



-)2 



the dying-out behavior of the resonance helicity amplitudes in the Q -dependence. Also, it was 
shown that a smaller value a = 331.7 MeV with Fg^Q"^) = 1 could nicely reproduce the helicity 
amplitudes for the 5*11(1535). See the dot-dashed curve in Fig. Ic of Ref. 
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explicitly included, i.e., 5ii(1535), L)i3(1520), Pi3(1720), Fi5(1680), Pii(1440), Fii(1710), 
Pi3(1900), and Fi5(2000). In the SU(6)(S)0(3) symmetry limit, the relative strengths of 
all these resonances and higher mass terms are fixed with the same quark-meson coupling 
denoted by a,,, apart from the quark model couplings. The configuration mixing among 
excited states, which breaks the SU(6)®0(3) symmetry, in principle should be taken into 
account. However, since the complexity arising from such a scheme somehow cannot be 
easily controlled, one can take an empirical strategy to introduce the strength coefficient 
Cn* for each resonance [p|,p3|. These coefficients can be determined by numerical fits, and 
their deviations from unit will refiect the breaking of the SU(6)(8>0(3) symmetry. In this 
study, the selected kinematical region allows us to focus on the two resonances S'ii(1535) and 
Di3(1520). Thus, only Cs^^ and Cd^^ introduced and treated as parameters. It should be 
noted that, coefficients Cs^^ and were found close to unit in a recent photoproduction 
study . 



Finally, we include an additional coefficient Cg^^ for the the longitudinal excitation am- 
plitude of the 5*11 in the electroproduction process. This quantity will be discussed in the 
following Section with respect to the experimental indications P5|,PB|. 



B. Analysis of observables 

With those parameters (Table D extracted in the data fits, we shall proceed to the 
detailed analyses of the model-fitting results. Predictions thus can be made for various 
observables. Some dynamical features arising from the electroproduction processes will be 
investigated. The 5*11(1535) properties can be reasonably constrained. 



1. Cross- sections 

In Fig. 0, the fitting results for the photoproduction process are shown for energies from 
threshold up to E^^^ = 900.4 MeV. The full curves come from the model outlined above and 
reproduce well enough the data. The dashed curve in Fig. I at P^ = 790 MeV denotes the 
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result in the absence of the 1^13(1520) contribution at the 5*11(1535) energy. Comparison 
between the full and dashed curve shows that the Di3(1520) plays a small but non-negligible 
role in 7p — > rjp. It accounts for the main non-S-wave behavior in the angular distributions. 
Meanwhile, the Su dominance is displayed by the (almost) isotropic behavior of the dashed 
curve. 

In Fig. ^ the differential cross sections for the electroproduction process at low-Q^ are 
presented (full curve) and compared with the old data at =0.4 (GeV/c)^ with e = 0.79 
and 0* = 0° ||2^. Components of the cross section (Eq. |l^) are also presented: transverse 
(dashed curve), longitudinal (dotted curve), transverse-transverse (TT) interfering (dash- 
dotted), and transverse-longitudinal (TL) interfering component (heavy dotted curve). Since 
the longitudinal contribution is quite small, it is the TT interference that accounts for 
the structure at small angles as shown by the solid curve. Meanwhile, it shows that the 
data prefer even smaller longitudinal cross sections at this region, although the fits are in 
agreement with the data reasonably. 

The recent precise measurements [^] of the t] meson electroproduction at =2.4 and 
3.6 (GeV/c)^ are successfully fitted by this approach. In Fig. ^, the angular distributions 
at W = 1.54 GeV with different azimuthal angles 0* for both momentum transfers are 
presented. The polarization parameters are e = 0.51 for = 2.4 (GeV/c)^ and 0.46 for 

= 3.6 (GeV/c)^. The solid curves represent the results for = 3.6 (GeV/c)^ and 
the dashed curves for = 2.4 (GeV/c)^. The results in the absence of the i5i3(1520) 
at = 3.6 (GeV/c)^ are shown by the dotted curves. It confirms that the non-S'-wave 
feature of the data is due to the Di3(1520) interference. Meanwhile, other resonances and 
the nucleon pole terms have only negligible influences. At = 3.6 (GeV/c)^, the fits start 
to underestimate the experimental dottes. It might suggest that the NRCQM form factor 
for the 1^13(1520) becomes inappropriate with the increasing Q^. In Ref. |^^, the data 
showed that the D13 form factor had more steeper slope with the increasing Q^, which led 
to a negligible D13 contribution above = 3.0 (GeV/c)^. However, in Fig. |, the D13 still 
makes sense at = 3.6 (GeV/c)^, which means it does not decrease as fast as the data 
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would require. 

With the parameters fixed, in Fig. the Q^- dependence of the transverse and longi- 
tudinal components of the total cross section are shown by the dashed and dotted curves, 
respectively. We plot the total cross section atot = <7t + (full curve), at the W = 1.535 
GeV, with £=0.6. The value of the latter quantity corresponds to the JLab experiment 
kinematics ||15|. The data are correctly reproduced. Although the values of e, determined 
by the experimental kinematics, varies from one datum set to another, the chosen value for 
e is not crucial in the illustration of the calculation. This is due to the fact that ai comes 
out much smaller than o"t. 

An interesting quantity to be investigated is the ratio RiQ"^) = ai/ax- Its Q^-dependence 
was found to be quite smooth in both experiment [p5| , p6| , p8| and theory . 



As a preparation, the two well known relations for an excited resonance in 7(i,)iV rjN 
are presented: 

and 

where, T ^ is the total width of the resonance and hr, = Tn/Tji is the branching ratio of the 
resonance decaying into rjN . 

In the calculations, we find that the inclusive (all contributions) and exclusive (only Su) 
cross sections produce the same ratios because the Su dominates in both ai and ar- Then, 
the ratio can be expressed as 

«"(«^) - i " WW? ■ 

2 

Note that, the amplitude Az vanishes in the Su excitations. 

2 

However, in both inclusive and exclusive calculations, the present NRCQM study slightly 
overestimates the experimental result Empirically, we define a parameter Cg^^ for the 



^ii longitudinal amplitude such that 
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Re.p{Q') = {C^J'RthiQ') ■ (33) 



Given R^xp = 0.23±0.11 at = 0.4 {GeV/cf from Ref. |3§, we find C^^ = (Rexp/Rth)^ = 
0.70. This result is in good agreement with the fitted value C^^^ = 0.65 in Table |, and 
explains the introduction of this parameter for the 5*11. The successful extraction of Cg^^ in 
the fits suggests that the leading order relativistic correction to the 5*1 could be absorbed 
into the parameter Cg^^ and such a correction might not be sensitive to a wide range. 
The latter point can be explained by the flat dependence of the ratio. In Fig. ^ the 
ratio RthiQ^) = ctl/o't is calculated using different values for Cg^^. The sohd curve denotes 
the result for = 0.65, the dashed for = 1, and the dotted for = 0.70. The 

i>il ' ill ' bii 



maximum ratio is at ^ 0.5 (GeV/c)^, which is in agreement with the finding of Ref. p8 



In the following part, we will see that the overestimation of the longitudinal cross section 
ai comes mainly from the photon vertex, where the longitudinal helicity amplitude Si has 
been overestimated. 



2. Sii(1535) photo- and electro- excitation helicity amplitudes 



To better understand the overestimation of the ai, we independently calculate 5*1 of the 



transition 7(i,)p — ^ 5*11(1535) in the EVF, 



^7 



(34) 



where l/jk is the Lorentz boost factor defined in previous Section. 

On the other hand, we explicitly write out the longitudinal cross section crL{Su) for the 
^ii(1535). 
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2 g-(kV7^W/7?)/3a2 
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(35) 



where the Lorentz boost factors have been included and the kinematical condition s = 

at the 5*11(1535) cm. energy has been used. Substituting the above equation into Eq. |3l|, 

5*1 can be also derived: 

2 



15 



^ 2 ^ aea5ii|q||k| 2 



2 



1 
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(36) 



Ef + Mat 3a27g_ 

In Fig. ^, the results for 5*1 from Eq. ^ (dotted curve) and from Eq. ^with Cg^^ = 0.65 
(solid curve) and C^^^ = 1 (dashed curve), are presented. Interestingly, it shows that 
although the dotted and dashed curves both overestimate the helicity amplitude Si, their 
values are very close. We thus equate these two expressions to give: 

bn = ^X, (37) 



Su 



where 



X = 



2 -"ii^ 



e 
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(38) 



Ml 'E/ + M^ '3a2^g_ 

and as^^ = a,,(Cs'iJ^ is the rjpSu couping constant. 

Intuitively, the above expression is a reasonable deduction since Eq. |3^ is exactly the 
branching ratio obtained by calculating 5*11(1535) rjp. However, recalling that fe^ and 
r^^^ are quantities to be investigated here, relations (as Eqs. |T7| - |55D derived from only 



5*11(1535) TIP cannot help us much, not mention that 0:5^^ is also unclear. We thus expect 
that information about these quantities can be derived in the photo- and electroproduction 
process. Unfortunately, another difficulty arises from the electromagnetic vertex in photo- 
and electroproduction. Since the experiment cannot separate the '-fNSu from the strong 
rjNSu couplong, studies of Su and 5*11 rjN are strongly model-dependent pT 



Therefore, the fit of the photo- and electroproduction data generally tells us what is the 
"parameter" (e.g. as^i) for 5*11 rjN instead of what is the coupling for 5*11 rjN . 

The overestimation of the 5i is an example exposing this kind of problems. 

Now, the trivial-look Eqs. |3^ and ^ can tell us more about the rjNSu coupling. First, 
since these quantities, 6^, Fg^j, and a^j^, are determined by fitting the reaction, we thus 
should bear in mind that their values have contained uncertainties from the electromagnetic 
interactions. Also, the model-dependent feature arising from non-5'ii contributions could 
mix into the fitted values for those parameters. As shown by Fig. ^ although the fitted 
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quantity Cg_^_^ = 0.65 better accounts for the data, we obviously cannot equate the sohd 



curve to the dotted one. In that case, Equations |3^ and |38| would not be a natural deduction 
at all. The non-trivial thing occurs only if we take all parameters fitted by the data except 
for Cg^j = 1 to produce the dashed curve in Fig. |^. Namely, both overestimations of the 
electromagnetic coupling turn to be equal to each other. Thus, it leads to the relation of 
Eq. in which the electromagnetic coupling has been cleanly removed. Since on the one 



hand, no information about the electromagnetic interaction is contained in Eqs. ^ and p8 
any more, and on the other hand, all quantities involved are determined by the fits, these two 
equations in effect provide an additional relation for all those fitted quantities. We thus can 
ask the question concerning the self-consistence of the model fitting results, whether those 
quantities fitted by the reaction data can satisfy the relation^ Again, we have to bear in mind 
that those quantities could have contained all uncertainties arising from the electromagnetic 
coupling as well as model-dependent features from the non-S*!! terms. 

To answer the question, we calculate the partial decay width of the 5*11(1535) using those 
fitted quantities, i.e. = ^r^Tsn, and F^ = asi^X; where x is calculated by Eq. The 
results are presented in Table |T[ The PDG estimations are also presented. It shows 



that F^ and F^ are in good agreement with each other, although discrepancies also exist. 
Quite clearly, the discrepancies reveal the uncertainties arising from the model-dependent 
aspects in the fits, while the consistencies highlight a reasonable constraint to the strong 
coupling. It suggests that the effective Lagrangian successfully accounts for the rjNSu 
coupling within an accuracy of 15%. We shall see below that this conclusion is important 
for the assessment of the 5*11 transverse helicity amplitude derived in the reaction. Note 
that the harmonic oscillator strength a has a moderate value 384.5 MeV, which is within 
the range of commonly used values, 330-430 MeV, this provides another justification for this 
approach. In this sense, the relation of Eq. |3^ could be regarded as being "satisfied" rather 
than "fitted". 

In the following part, we turn to the results for the transverse helicity amplitude for the 



S'ii(1535). From Eq. the transverse hehcity amplitude for the S'ii(1535) can be expressed 
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as 



2Mn brj 

n aeOsiilql {Ef + Mn] 



2M, 



X + 7^ 



N 

2 g-(kV7^+qV7,2)/3a2 



LOr. 



+ 1 



2q2 



(39) 



9 [nq \Ef + MN ' 7 3a27^_ 
where crT(7(i.)P — > rjp) is the exclusive transverse cross section calculated by the model, 

and the Lorentz boost factors are taken into account. The result of Eq. ^ is presented 
in Fig. 1^ (dashed curve), which is in good agreement with the experimental data. It is 
worth noting that a direct relation as Eq. |^ cannot be obtained for the transverse helicity 
amplitude. The basic reason is that the EM operators used in this model [0 are essentially 
different from the NRCQM ones However, as shown by the dashed curve in Fig. |^, 

those well constrained parameters indeed provide a reasonable discription of the S'ii(1535). 



Such a result certainly highlights again the relation provided by Eq. ^ and It is also 
useful for the understanding of a variable defined in the literature. In Ref. 0, it was found 



that the quantity = y^^riAi/Tsu, where is the phase space factor, was not sensitive 
to the fitting scheme since uncertainties within and F^n went to the same direction. 
That feature arises because the constraint to the riNSu vertex was not enough. However, 
here we find that Eq. ^ restricts the 6^ and Fg^^ in an inverse direction compared to C,t- 
Supposing and Fg^^ both increase or decrease, the relation of Eq. will be destructively 



unbalanced. The satisfactory of Eq. |33 thus indeed serves as a further constraint to those 
fitted quantities. 

The dominant S'ii(1535) production can be investigated by substituting the exclusive 
(^t{1{v)P Sii rjp) with the inclusive a-r in Eq. 



Ms,, Ts,, 
2Mn br, 



ax ■ 



(40) 



The difference between Eqs. p9^ and 40 reflects the influence of the non-5'-wave contributions 



in the transverse helicity amplitude. The calculations of Eq. |3^ (dashed curve), and Eq. ^ 
(solid curve) are presented in Fig. ^. It shows that the impact of the non-^n contributions 
produces about 3% difference in the helicity amplitude, and the transverse cross sections are 
dominated by the S'ii(1535) excitation. 

In Ref. the impact from the longitudinal cross section have been investigated through 
the following relation: 



12 c^t(7P Sn r]p) 



' ^' 2M;v b, (l + £i?(Q2)) • ^ ^ 

As shown in Fig. ^ the ratio -R(Q^) is much smaller than unit, thus, the effect of the 
longitudinal contribution, reflected by the factor 1/(1 + eR^Q"^)), is found negligible in 
comparison with Eq. ^ especially at high regions. 

In Fig. ^, we also present the calculation of Ai based on the dipole-like Lorentz boost 
factor in the meson-nucleon cm. frame (Eq. for the S'ii(1535) (See the dot-dashed 
curve). With the same set of parameters, the dot-dashed curve apparently under-estimates 
the experimental values, although compared to a parametrized dipole form factor, 1/(1 + 
(5^/0.8)^ (see the dotted curve), it better accounts for the trend of data. 

The resonance contributions to the transverse helicity amplitude Ai(Q^) drop down 
with the increasing as shown by the full curve in Fig. |[ At high region, the pQCD 
dominance in the helicity-conserved amplitude (Q^) predicted the asymptotic behavior of 
Ai (Q^) fall-off like l/Q^ |3^ . That is to say, the quantity Q^Ai (Q^) will approach a constant 
at high regions. However, the studies by Isgur and Smith ||4^, and Radyushkin |Q led 
to quite different results. They found that such a scaling behavior in the exclusive processes 
were still dominated by the non-perturbative contributions rather than the perturbative 
ones over wide regions. We present the calculation of the quantity Q^AilQ"^) in Fig. ^to 
compare with the data. At = 2.4 and 3.6 (GeV/c)^, the quantity Q^Ai{Q'^) is calculated 
through Eq. It gives 0.201 and 0.268 GeV^, respectively, which are consistent with 
the non-scaling behavior found in experiment at small Q^. At Q"^ = 4.0 (GeV/c)^, we 
get Q^Ai = 0.282 GeV^. This quantity approaches the maximum 0.302 GeV2 at = 
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5.5 (GeV/c)^, and then falls down slowly with the increasing . However, due to the 
shortcoming of the NRCQM, the results above 4 (GeV/C)^ should not be taken seriously. 
More fundamental approach is needed for the dying-out region of non-purterbative processes, 
where a continuing scaling behavior from pQCD processes might emerge. 



IV. SUMMARY AND CONCLUSION 

This work was motivated by the recent electroproduction data from JLab [T^ in the 
dominant kinematical region. We also took advantage of even more recent photoproduction 
data from GRAAL 0]. Such a database, expected to be significantly enlarged in the near 
future, offers an excellent opportunity to investigate the reaction mechanism of the t] meson 
photo- and electroproduction. In particular, the selected kinematical region allows us to 
study the nature of the 5*11(1535) resonance with unprecedently precise data. 

The model presented here is an extension of a constituent quark model approach with 
a chiral effective Lagrangian [jl2[ to the electroproduction process. On the basis of the 
symmetric NRCQM, there are only a limited number of parameters appearing in this ap- 
proach. In principle, the breaking of this symmetry can be phenomenologically introduced 
through an additional coefficient for each resonance, which then could be determined in 



the numerical study [g,^]. However, different from Refs. [p|,p3|, in this work attentions 
are paid to the 511(1535) and Di3(1520). The selected kinematics permit us to introduce 
the symmetry breaking coefficients only for the 5ii(1535) and Di3(1520) due to the domi- 
nance of the former and significant interferences from the latter. For other resonances, their 
symmetry breaking coefficients are neglected due to their negligible effects. In this work, 
the 5*11 (1535) 's total width, branching ratio of riN channel, as well as the harmonic oscilla- 
tor strength of the quark model are treated as parameters and determined by the fits. The 
quark model symmetry breaking could introduce contributions from the resonance 5ii(1650) 
and Di3(1600), which have been suppressed in the SU(6)cg>0(3) symmetry limit due to the 



Moorhouse selection rule |T3] . Although an explicit consideration of such a breaking has not 
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been doen here, a possible influence of these resonances could have been contained in those 
fitted quantities as a background contribution. In this sense, the 15% uncertainty for the 
r/A^ 5*11 (1535) coupling can be regarded as a reasonable estimation. Another feature of this 
approach is that the Lorentz boost factor in the equal velocity frame (EVF) succeeds in the 
description of the 5*11(1535) form factor. Therefore, an ad hoc form factor is avoided. 

We have shown that the extracted parameters are compatible with their values coming 
from previous independent works. The cross section data for both photo- and electropro- 
duction are well reproduced. Explicitly, it shows that the electroproduction process indeed 
provides us with rich information. It sheds light on the r7A^S'ii(1535) coupling, which had 
not been well-determined in the real photon reaction. Although the scalar coupling of the 
longitudinal photon is found to be larger than that reported in Ref. ||35|, we succeed in 



removing the model-dependent uncertainties arising at the electromagnetic interaction ver- 
tex {■y{v)P Sii) from the transition amplitude 7(t,)P — > rip for the 5*11(1535). The 
derived simple relation (Eq. |37| ) thus provides a constraint on the rjNSu coupling with an 
uncertainty of about 15%. This could be the first explicitly constrained estimation of the 
5*11(1535) properties in theory. Such a result might lead us to the conclusion that the main 
component in the 5*11(1535) wave function is rather like a non-exotic three-quark state than 
a KA (KTi) bound state Possible configurations from the latter should be small at the 



5ii(1535) energy, if exists. 

It should be noted that this calculation was done before the publication of a new set of 
data from the CLAS Collaboration [|^. The 5*ii helicity amplitude Ai for 0.25 < < 



1.5 (GeV/c)^ has been significantly improved by that experiment. Their results, which are 
highly model-selective, shows in good agreement with our calculations. 
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TABLE I. 
the text. 



TABLES 

Extracted values for parameters of this model. Detailed discussions are given in 



a (MeV) Ts,, (MeV) 




Cdi3 




0.10 ±0.01 384.5 ±0.2 143.3 ± 0.2 0.55 ± 0.01 


1.20 ±0.30 


-0.92 ±0.03 


0.65 ±0.01 2.36 



TABLE II. Extracted total width and partial decay width T" = bf^Fsn and = asnX 
for the /S'ii(1535). Corresponding estimations by the Particle Data Group are also listed. 

r5n (MeV) n (MeV) (MeV) 



This work 143.3 ± 0.2 78.8 89.4 

PDG 150 75.0 
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FIGURES 
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FIG. 1. The Q'^-dependence of the Lorentz boost factor in the EVF (solid) and meson-nucleon 
cm. frame (dashed), respectively. 
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FIG. 2. Differential cross sections in the real photon limit. Solid curves denote the fitting 
results, while the dashed curve denotes the result in the absence of the Di3(1520). Data are from 
Mainz Q (triangles) and GRAAL [|] (circles). 
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FIG. 3. Differential cross section at = o.4 (GeV/c)^ with (p*=0°. The solid curve comes 
from the global fits. Different components of the cross section are also presented: the dashed curve 
for daT/dQ*; the dotted for daL/dQ*; the dot-dashed for daTr/dfl*; and the heavy-dotted for 
daTh/ d^* ■ Data are from Ref. p9| . 
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FIG. 4. Fitting results for the differential cross sections at = 2.4 (dashed) and 3.6 (GeV/c)^ 
(solid) for different (/>*. The dotted curves are calculations without Di3(1520) contributions at 
q2 ^ 3 g (GeV/c)2. Data are from Ref. [|T5[. The triangles are data at = 3.5 (GeV/c)^ while 
the full dottes are data at = 2.4 (GeV/c)^. 
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FIG. 5. The Q'^-dependence of the total cross sections (sohd curve). The total transverse and 
longitudinal cross sections are also shown by the dashed and dotted curves, respectively. Data 
come from Refs. [|J2|M|l5lEnj3|] . 
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FIG. 6. The Q^-dependence of the ratio of a^/aT- The solid and dotted curves correspond to 
Cg^^ = 0.65 and 0.70, respectively, while the dashed curve is for Cg^^ = 1. The full circle is from 
Ref. |^5| and the triangle from Ref. ||36| . 
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FIG. 7. Longitudinal helicity amplitude Si for S'ii(1535). The solid curve and dashed curves 

2 

are calculated with C^^^ = 0.65 and Cg_^^ = 1, respectively, by Eq. while the dotted curve is 



given by Eq. 34 
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FIG. 8. Transverse helicity amplitude Ai for the S'ii(1535). The solid and dashed curves 

2 

are calculated by Eqs. ^ and respectively. The dot-dashed curve is calculated based on the 
Lorentz boost factor in the meson-nucleon cm. frame, and the dotted curve based on a dipole 
form factor. The full circles are from Ref. |15] and other data points from Refs. |l,37,|^. 
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FIG. 9. The quantity AiiQ"^) calculated in this model. The full circles are from Ref. p!5 

2 ^ 

and the empty ones from Refs. [p],|37|,pt . 
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